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We develop a scaling ansatz for the master equation in Dvali, Gabadadze, Porrati cosmologies, 
which allows us to solve the equations of motion for perturbations off the brane during periods when 
the on-brane evolution is scale-free. This allows us to understand the behavior of the gravitational 
potentials outside the horizon at high redshifts and close to the horizon today. We confirm that the 
results of Koyama and Maartens are valid at scales relevant for observations such as galaxy-ISW 
correlation. At larger scales, there is an additional suppression of the potential which reduces the 
growth rate even further and would strengthen the ISW effect. 



I. INTRODUCTION 

That cosmic acceleration is a fact appears indubitable. 
Instead of an exotic new form of dark energy driving the 
acceleration, it may be caused by a modification of grav- 
ity. Precise measurements for gravity are only available 
in the range of scales from a millimeter to that of the so- 
lar system — we do not have any direct probe of Einstein 
gravity beyond these boundaries. Cosmic acceleration 
may originate in a breakdown of Einstein gravity at dis- 
tances beyond the range above. 

Dvali, Gabadadze and Porrati (DGP) [Hi have pro- 
posed a braneworld theory in which our universe is 
a (3+l)-dimensional brane embedded in an infinite 
Minkowski bulk. Gravity propagates everywhere, but, on 
the brane, an additional (3+l)-dimensional gravitational 
interaction is induced. This allows for gravitational po- 
tentials on the brane of a (3+l)-dimensional form at 
small distances to evolve into (4+l)-dimensional form be- 
yond a crossover scale determined by the unknown energy 
scale for the bulk gravity. The cosmological solution of 
this theory was shown to exhibit accelerated cosmic ex- 
pansion without the aid of an exotic energy component 
like dark energy 01 0. 

It has already been shown that the linearized field the- 
ory as defined by the DGP model contains ghost degrees 
of freedom [^] [5J [Sj [Ij [SJ , or even may violate causality in 
certain limits [9||. It is known, for instance, that the de 
Sitter background is unstable to classical linear perturba- 
tions; however, it is claimed in that strong-coupling 
effects at small radii around matter sources ensure that 
the theory remains stable. 

The point of view of our work is to assume that linear 
perturbation theory remains valid on the largest scales. 
This is motivated by the fact that the late universe 
is dominated by the gravitational interaction of dark- 
matter haloes. The internal structure of the haloes is 
controlled by the strongly-coupled non-linear theory. On 
the other hand, the radius below which strong-coupling is 
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important for haloes is approximately equivalent to their 
size and therefore their interactions should be driven by 
the linear theory analyzed in this work. 

We do find that deep into the accelerated era the space- 
time becomes unstable on the timescale of the expansion. 
However, this is an efifect that only becomes important 
far into the future and is negligible as far as the observa- 
tional impact today is concerned. We therefore assume 
that during the early universe, when the theory does not 
exhibit instabilities, the analysis for DGP proceeds in ex- 
actly the same way as that for GR. Then, deep during 
the acceleration era, instabilities develop and the theory 
may or may not be saved by non-linear eflects — an issue 
on which we remain agnostic. This evolutionary history 
appears to be the only one which is capable of reproduc- 
ing the universe as we see it. If strong coupling effects 
are important straight away and at all scales, the ap- 
proximation of a homogeneous background cosmology is 
completely inapplicable and the DGP model would not 
be able to reproduce observations such as supernova lu- 
minosities. We therefore effectively assume a best-case 
scenario for DGP: should this analysis fail to predict the 
observations, the model is excluded. If it passes the ob- 
servational tests, a more careful study of the effects of the 
strong coupling regime during the acceleration era would 
be required to understand fully the future evolution of 
the universe in the DGP model. 

Under the above assumptions, the equations of mo- 
tion for the theory of gravity on the brane, pertinent to 
the study of cosmology, do not close owing to the inter- 
action with the bulk at first order in perturbation the- 
ory. Koyama and Maartens have used a quasi-static 
approximation, valid well within the horizon, to investi- 
gate structure formation at smaller scales. This solution 
shows the essential role that the bulk plays in correcting 
the gravitational potentials, reducing the growth rate. 
Also, Lue et al. have reached a similar conclusion 
using a diflerent approach, including non-linearities in 
their calculations. 

In the following, we present a new scaling ansatz for 
the master equation, allowing us to solve the equation 
and calculate the resulting cosmological evolution at all 
scales for high redshifts, and close to the horizon today. 
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In we review the linearized equations of motion for 
DGP. We present our scaling solution in §1111 and dis- 
cuss its cosmological implications in §IVI We study the 
robustness of the scaling solution in the Appendix and 
discuss these results in ^Vl 



II. DGP EQUATIONS OF MOTION 

A. Background 

In the DGP model, gravity alone propagates in the 
bulk, and the 5D gravitational theory is complemented 
by an induced 4D Ricci scalar restricted to the brane. 
We assume that both the bulk and the brane have zero 
tensions, i.e. the cosmological constants are zero. We 
thus start off with the basic DGP action: 
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a 3-brane embedded in an empty bulk, with all the 
standard-model fields localized on the brane at x = 0. 
The constants fi^ and define the energy scales of the 
theories of gravity: one is Newton's constant, fi^ = 87tG, 
the other represents the energy scale of the bulk gravity. 

As shown by Dvali et al. the ratio of the two 

scales defines a cross-over radius beyond which the four- 
dimensional gravitational theory transitions into a five- 
dimensional regime 
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This scale is chosen to be of the order of the current Hub- 
ble length so that the acceleration of the expansion to- 
day results from the 4D to 5D transition. For illustrative 
purposes, we will take in our calculations TcHq = 1.32, 
in a matter and radiation universe with fim — 0.24 
and h = 0.66. When comparing with ACDM, the con- 
cordance set of parameters will be used: flm = 0.25, 
h = 0.72. The two sets of parameters represent cos- 
mologies which are best fits to supernova luminosity data 
(SNLS and the distance to the last-scattering sur- 
face (WMAP but assuming that the universe is fiat. 
We discuss the fits to these data in a companion paper. 
Song et al. 

For the purpose of cosmological calculations, we would 
like to reduce the five-dimensional braneworld to an ef- 
fective theory on the brane, which could be studied us- 
ing the usual range of four-dimensional tools. Using a 
4+1 decomposition of the theory [lB| [l3| [3 j can de- 
rive the effective on-brane equations of motion — a set of 
modified Einstein equations 



(3) 



where the Greek indices range across all the dimensions, 
= 0,1,..., 4. f^i, is a tensor quadratic in the 4- 



dimensional Einstein and energy-momentum tensors, 
1 1 if A^ 

fliu = -y^AA^j^i, — -A^^A^a + -^9tiu [AafiA'^ — 
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while Efj^^ is the bulk Weyl tensor projected onto the 
brane using the vector normal to the brane, 



(6) 



It can be shown that for branes with maximally- 
symmetric spatial hypersurfaces, the projected Weyl ten- 
sor must take the form Ca~^, with C a constant (see [l9j : 
it is effectively a constant of integration for the back- 
ground). Since the Weyl tensor is traceless, it will dilute 
as rapidly as radiation and will become irrelevant at late 
times, if it doesn't completely dominate the dynamics 
initially. We will therefore set C — 0, allowing us to find 
the modified Friedman equation for the background evo- 
lution of the cosmology on the brane. In the case of the 
fiat brane, which is the only one to be considered here. 
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leading to the result that, for the upper-sign selection, 
the cosmology tends to a de-Sitter phase as the matter 
density gets diluted away, potentially providing a model 
for the observed acceleration [23|. This is the choice we 
will make henceforth. 

In addition to Tc, there is another scale present in 
the theory, which thus far we have ignored: the strong- 
coupling scale given by 



1/3 



(8) 



where Tg is the Schwarzschild radius of the point source 
under consideration. Beneath this radius, the linear ap- 
proximation developed above is not valid and the theory 
returns to 4D GR with small corrections. It is not un- 
derstood how two sources superimpose when the linear 
regime is not vaHd. However, the Hnear density fields 
are actually constructed out of a spatial average of dark- 
matter halos, each of which has r* comparable to its ra- 
dius. Thus we should be able to use the linear theory to 
discuss the gravitational dynamics of the spatial averages 
of the halos, but most Hkely not their internal structure. 



B. Linear Perturbations 

An observer restricted to the brane will perceive the 
universe as purely four dimensional. Therefore, the most 
general linear scalar perturbations to the induced fiat 
four-dimensional metric can be parameterized by 
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while he linearized energy-momentum tensor can be writ- 
ten down as 



T° = (1 + w)pd^q , 



(10) 

(11) 



Here p is the density of the cosmological background, 
w = p/p is the background equation of state parame- 
ter, and Cg = Sp/Sp is the sound speed for the pressure 
perturbations. 

The full 5D perturbations of the bulk are richer: Def- 
fayet in [23| has shown that their effect on the brane 
can be reduced to the presence of perturbations to the 
Weyl tensor as an additional source of stress-energy. The 
precise relationship between the various components of 
the Weyl tensor is determined by their relationship to 
a master variable (see ^11 Cp . which in turn encodes all 
the gauge-invariant perturbations to the full 5D metric 
throughout the bulk. We define the linear scalar pertur- 
bations to the Weyl tensor by 



E% = -fi^pSE , 

E\ = p^pd.QE 
A ■ = ^ WEP 



SeS% 



d'd, - -5\d' ] ^E 



(13) 
(14) 

(15) 



We have here assumed adiabatic pressure fluctuations in 
the multicomponent matter system Cg = p/ p. 

On the other hand, the Weyl tensor is not separately 
conserved and its equations of motion come from the 
Bianchi identity, V^G^^ = 0, as appHed to Eq. ([3]), 



(12) V-i?^, = 4r2v^^, 



We can rewrite this as 

SE + {l~iw)HSE-^qE^O, 

1 2 , 
q'e ~ SwHqE + -Se - gfc TTE = S , 

where the source term 
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Thus a non-zero Weyl tensor is unavoidably generated by 
matter perturbations in Hnear theory [ll,]. 

In order to close the above equations, we need the 
analogue of an equation of state to relate the Weyl 
anisotropic stress tt^; to the other components of the Weyl 
tensor 5e and qE- Unlike the relationship between the 
Weyl pressure and energy density, this relation requires 
a consideration of perturbations in the bulk. 



The Weyl tensor is traceless, and hence the pressure per- 
turbation will behave like radiation, i.e. we — 1/3. The 
linearized Poisson and anisotropy equations then are 
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with A the comoving density contrast and the equivalent 
definition for Ae'- 
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The equations of motion for the energy momentum 
tensor are given by the conservation law V^T^^ = sup- 
plemented by equations of state that define the stress 
fiuctuations 
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C. Master Equation 

In [2l|, Mukohyama showed that for maximally- 
symmetric five-dimensional space-times, the full five- 
dimensional linear scalar perturbations in the bulk can be 
described using a master variable f2. Since the bulk being 
considered is just Minkowski, and the brane is assumed 
fiat, we can parameterize the unperturbed background 
5D metric by [H 



where the brane sits at y = and 



a{l + H\y\), n 
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This parameterization is also valid for branes with non- 
zero curvature provided that \ f2k \ 1- The master vari- 
able then obeys a hyperbolic equation of motion: 




„ , n df2\ nk'^ 
TT- - -T^fi = 0. 
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We can then express all the gauge-invariant perturba- 
tions to the 5D bulk as functions of derivatives of the 
master variable. In particular, Deffayet has shown 
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that the components of the Weyl tensor evaluated on the 
brane can be expressed as 




(28) 
(29) 
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We will hereafter implicitly assume evaluation at ?; = 
for the master variable in the on-brane equations where 
no confusion might arise. We can now rewrite the Bianchi 
identity in terms of the master variable, obtaining, after 
assuming that the cosmological fluid has no anisotropy 
stress, 



n-3HF{H)f2 
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where R expresses the derivative across the brane 
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We have defined two new functions of the Hubble param- 
eter 
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Were it not for the dQ/dy derivative across the brane in 
i?, this equation would be a simple dynamic equation for 
J?, which, given an evolution equation for the source Z\, 
could be solved as a coupled equation. The role of the 
master equation is to define R, the relationship between 
dfljdy and fl. 

Koyama and Maartens [lH adopted a quasi-static ap- 
proach to solve these equations. The master equation 
then implies that the gradient 
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In the Bianchi identity the time-derivative and brane- 
derivative terms are neglected compared to those of or- 
der (k/aHYQ. This quasi-static approach leads to their 
solution to which we will refer henceforth as "QS" 
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This solution is equivalent to a closure relationship for tte 
in terms of 6e through Eqs. (|28l) and (|30| . Using this, 
we can define the quasi-static limit (QS) of DGP gravity 
where the Poisson and anisotropy equations become 
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These equations are equivalent to the linear limit of the 
results obtained by Lue et al. [l^ . 

In the next two sections, we will show how the quasi- 
static solution is dynamically achieved for the perturba- 
tions shortly after horizon crossing and discuss large-scale 
deviations from this solution. 



III. SCALING SOLUTION TO MASTER 
EQUATION 

A. Causal Horizon 

The master equation is a wave equation sourced by the 
comoving density perturbations on the brane through the 
Bianchi identity (|3T|) . 

Since, in appropriate coordinates, the bulk is just 
Minkowski, the evolution of 12 in the bulk can be seen 
as a normal propagating wave given a boundary condi- 
tion from the behavior on the brane. Beyond the causal 
horizon, the bulk should remain unperturbed. This 
causal horizon must be invariant in all coordinizations 
of the bulk; therefore we can locate it by finding the null 
geodesic of Eq. l(25|) . giving us the y-position of the hori- 
zon as 



i = yhorH = aH 



da 
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Before the acceleration epoch, this reduces to ^ = 1/(2 + 
3w) for a cosmology with a constant equation of state pa- 
rameter, i.e. ^ = 1/3 for a radiation-only cosmology and 
^ = 1/2 for a matter-only cosmology. We are making 
the assumption that the universe has not gone through a 
period of infiation, which would have moved the horizon 
much further out. If infiation did take place, it is not un- 
reasonable to expect that any perturbations that existed 
prior to the infiationary phase will have been pushed far 
away and the bulk will start in an unperturbed state in 
the vicinity of the brane at the beginning of radiation 
domination, resulting in a causal horizon equivalent to 
that of the cosmology with no infiation. 

The constancy of the horizon during the domination 
of a particular fiuid suggests that we can define a new 
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variable, in which the horizon will remain fixed at all 
times, lying at a; = 1 



yH 



(41) 



We can then recast the master equation in x and solve it 
as a boundary value problem with the value of J7 at the 
horizon, a; = 1, set to zero. 



B. Scaling Ansatz 

The second boundary condition needed to solve the 
master equation is the behavior of the master variable 



on the brane. During epochs when the source remains 
scale-free and the Bianchi dynamics also do not change, 
one would expect that the master variable also obeys 
a scaling ansatz on the brane f2\y=o — AdP , where A 
and p are constants. Likewise, during such epochs we 
expect the master variable in the bulk to reach a stable 
solution in the variable x, the distance in units of the 
causal horizon, for a given wavenumber k/aH . 

We therefore propose a new ansatz for the solution to 
the master equation (|27|) : 



n ^ A{p)aPG{x) . 
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With this assumption, the master equation l(27|) becomes 
the ordinary differential equation 
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where the primes denote differentiation with respect 
to In a — the new time coordinate which will be used 
henceforth — and h = H'/H. In deriving the above, we 
have neglected time derivatives of p and : the scaling 
ansatz is not expected to be valid when p is not a con- 
stant, i.e. during times when the cosmology is undergoing 
a change from the domination of one fluid to another. In 
addition, strictly speaking, G is actually a function of 
both X and k/aH, even in the scaling limit. We have 
therefore also assumed that k/aH is a constant, which is 
vaHd in the fc = Hmit. As we explain later, as k/aH 
approaches unity, where its time-derivatives might im- 
pact the solution significantly, the character of Eq. l(43|) 
changes and terms which do not involve the derivatives 
of k/aH dominate. 

Note that, provided w is constant, one of the denomi- 
nators in Eq. (|43l) can be re-expressed as 



l + x^{l + 2h) = l-x. 



(44) 



Thus the equation has a regular singular point x = 1, 
exactly at the causal horizon. This is not a coordinate 
singularity (all the entries of the metric are regular there) , 
but is a refiection of the junction between perturbed and 
unperturbed space-times. 

Supplying 77, p as a function of the scale factor is 
enough to solve this ODE as a boundary-value problem, 
requiring that G{x) be 1 on the brane and at a; = 1. 
This in turn gives the value of R as 
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and closes the evolution equations for the perturbations 
on the brane. 

We will henceforth refer to this solution as the dynam- 
ical scaHng or just scaling solution and use the acronym 
"DS". 



C. Iterative Solution 

In practice, one does not know the scaling index p a pri- 
ori and, moreover, it can change during the evolution of 
a k mode as the master variable leaves one scaling regime 
and enters another. We therefore solve for p iteratively 
by demanding consistency with the Bianchi identity. 

To determine the zeroth-order solution for p, we sub- 
stitute the ansatz Eq. (|42l) into Eq. ((3T|) to obtain 



A{p) = 



2a^-P 
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K{H)Hr, Hr, 
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Before the acceleration epoch, when the expansion is 
dominated by a single fiuid, J{H) and F{H) are constant, 
while K{H) is a simple power law in a. We therefore set 
p =p^°\ the zeroth order solution for k/aH <C 1 



P 



(°) = 3 + 



dln[pZ\(0)/(ii:(iJ)ij2)] 



(47) 



x=0 



d In a 

Note that this definition allows transitions between seal 
ing regimes where w changes. We shall see that the solu 
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tions for R become independent of p for k/aH ^ 1 and 
so we use this as the zeroth order solution for all modes, 
both super- and subhorizon. 

Finally we need a zeroth order solution for A^^^ . Modes 
of interest to large-scale cosmological tests are superhori- 
zon during radiation domination and enter the horizon 
either during matter domination or the current acceler- 
ation epoch. Before the acceleration epoch and in the 
absence of anisotropic stress, these modes obey [131 
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D{D + 1) 



(48) 



where D = a/aeq. In order to obtain the zeroth-order 
solution for p'-'^-' we assume that the growth prescribed 
by equation Eq. I|48p continues until today (see Fig.[T]). 

Given this, we solve the master equation Eq. (|43l) to 
obtain the off-brane gradient, R, and then dynamically 
solve the Bianchi identity Eq. (|3T|) coupled to the cosmol- 
ogy, for the particular mode. Once the dynamic evolu- 
tion for the first-order solution J7*^^-' is obtained, we can 
iteratively improve our estimation of p by numerically 
calculating 



dlnf^W 
din a 



(49) 



and repeating the above prescription. We find in the next 
section that this procedure converges quickly and alters 
the value of p only when p is not a constant, as expected. 
We display the effects of the iteration on p in Fig. [TJ 
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FIG. 1: Value of exponentp in the scaling solution f2\y^o « a,^ 
for superhorizon modes, plotted as a function of the scale fac- 
tor. The p'"' calculation assumes that the density perturba- 
tion A follows Eq. (gH) at all times. The p^*' result is the 
output of an iterative process, where p is used to calculate R 
and hence the evolution of fi which is in turn used to derive 
a correction to p. 



IV. COSMOLOGICAL IMPLICATIONS 
A. Limiting Cases and Numerical Solutions 

The evolution of the master variable exhibits sev- 
eral distinct phases that are distinguished by the on- 
brane scaling evolution: that during radiation domina- 
tion, matter domination and the de-Sitter acceleration 
phase both inside and outside the horizon. The scaling 
of the Bianchi identity Eq. ((3T|) determines the value of p 
at a particular scale factor, while the solution to the mas- 
ter equation Eq. (|43l) determines the off-brane gradient 
R. To better understand the nature of the solution and 
how it impacts perturbation evolution, we will compare 
the analytic expectations to the full numerical results in 
the various phases. 
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FIG. 2: Evolution of ratio of ofF-brane gradient to master 
variable, R, as defined in Eq. I|32p . for a selection of modes. 
On superhorizon scales, R is constant whenever p is a con- 
stant. Once the mode enters the horizon it rapidly approaches 
R = —k/aH. As the universe enters the de-Sitter phase, the 
modes again leave the horizon and R asymptotes to 1. 



1. Superhorizon Modes 

During radiation domination, Hvc ^ 1, H'/H ~ 
— 2 + a/2aeq and F{H) — 1/3. The particular com- 
bination in the denominator of K{H) causes the first- 
order contribution of H'/H to cancel, leaving us with 
K{H) oc a^^. The Bianchi identity Eq. |[3l]) for super- 
horizon modes then dictates that 



(50) 



Since A a: a? during radiation domination, this gives 
p = 6 for superhorizon modes. Inserted back into the 
master equation under the scaling ansatz Eq. (|43l) . this 
value of p implies 



R 



{k/aH <C 1, radiation domination) . (51) 
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With R determined, the equations of motion for the per- 
turbations on the brane are closed. 

This analytic expectation also serves as the initial con- 
ditions for the numerical scaling solution. In practice, we 
begin the integration at a = 10~®, when all modes of in- 
terest are outside the horizon. The numerical solution 
for p is shown in Fig.[T]and for R in Fig. [21 Note that, in 
the radiation-dominated era, their values stay stable at 
the analytic prediction for all iterations of the solution. 

The large-scale modes of interest remain outside the 
horizon during the whole radiation-dominated epoch. In 
the matter-dominated epoch, the evolution outside the 
horizon can be obtained by noting, F{H) = 1/2, K{H) — 
4, H'/H = -3/2, and Hrc > 1. The Bianchi identity 
dictates that 



A cx a^-PA . 



(52) 



Given that A (x a, p — A. Since the matter dominated 
solution is of particular interest, we explicitly give the 
master equation under the scaling ansatz 
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FIG. 3: Off-brane profile for G{x) obtained by solving equa- 
tion l|43p during matter domination (log a = —2), compared 
to off-brane profiles for the quasi-static (QS) solution. For 
high k/aH the profiles are very narrow and effectively inde- 
pendent of the position of the causal horizon: they penetrate 
very little into the bulk and the behavior of the solution is 
practically independent of the value of p. In this regime, the 
QS solution is practically coincident with the scaling solu- 
tion. For modes with low k/aH, the solution is non-zero in 
the whole interval x G [0, 1) and therefore it depends strongly 
on the value of p. QS severely underestimates the gradient of 
the profile in this regime. 



This is solved as a boundary value problem with bound- 
ary conditions G(0) = 1 and G(l) = 0. The form of the 
numerical solution to this equation is shown in Fig. [H In 
the large-scale Hmit, the gradient reaches 



R 



-1, 



{k/aH ^ 1, matter domination) 



(54) 



In the numerical solution of Figs. [T] and [21 these values 
of p = 4 and R = —1 are achieved gradually as the ex- 
pansion becomes matter dominated. The iteration of the 
numerical solution in fact further smooths the transition 
until a stable form is achieved as would be expected. Note 
also that R is very insensitive to k/aH, provided it be 
less than 1 such that the mode is larger than the horizon. 



2. Subhorizon Modes 

The modes of interest cross the horizon either during 
matter domination or the acceleration epoch. For large 
values of k/aH, the final term of the master equation 
Eq. (|43| dominates over other parts of the coefficients of 
G{x) and G'{x). As evidenced in Fig. [31 in this regime, 
the solution does not penetrate very far into the bulk. 
We can thus expand the master equation around x — 0, 
reducing it to 



d^G 
dx^ 



e(1^2^x)(^] G^O. 



(55) 



This matches the quasi-static approximation up to first 
order in x. Therefore, for k/aH ^ 1, the relation giving 
the gradient on the brane is exactly as in the QS ap- 
proximation [ll| , with no dependence on the value of p: 



R = — : {k/aH ^ 1, matter/acceleration). (56) 

aH 

The numerical solution for R evolved through horizon 
crossing is shown in Fig.[2l It reaches this scaling shortly 
after horizon crossing. 

Despite this independence of p, the master variable 
does achieve a scaling form during matter domination. 
The Bianchi identity l(3T|) can be reduced to the QS form, 
Eq. l(36l) . and implies 



Ai 



(57) 



and hence p = 3. 

In general then, i? is a function of both k/aH and p{a) 
and therefore each mode needs to be followed separately 
through its evolution both outside and inside the hori- 
zon. However, only outside the horizon does the value 
of the off-brane gradient actually affect the evolution on 
the brane, since for high k/aH it is subdominant in the 
Bianchi identity Eq. (|3T|) . 
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3. Asymptotic de-Sitter phase 

At late times, the DGP cosmology enters the self- 
accelerated de-Sitter phase. During this time all modes 
exit the horizon while the causal horizon in the bulk, 
^, grows rapidly toward infinity. This allows us to con- 
centrate on the k ^ aH limit. Deep into the de-Sitter 
phase of the expansion, Hrc = 1, K{H) = F{H) = 1. 
The master equation (|27|) can now be rewritten (in the 
scaling approximation) as 



dn p{p - 3) 



dy^ 1 + Hydy {I + Hyf 



i7 = 0. 



(58) 



This equation has an analytic solution. Assuming that 
the perturbations vanish at the causal horizon, J7 = at 
y = ^/i?, and that p ^ 3/2, we can find the off-brane 
gradient: 



«4(3-|2p-3|)- ,^Jg,;j,_^ 



(59) 



This can now be combined with the Bianchi identity 



(p2-3p+l + i?)i7 



(60) 



It can be shown that A becomes a constant during the 
acceleration era. 

In the limit where r^Z\/i7 — > and ^ — > oo, i.e. at 
very late times, the solution to the above two equations 
combined is i? = 1 and p = 2 or p = 1, where we would 
expect the fastest-growing mode, p = 2, to dominate. We 
plot the numerical solutions in Figs. [T] and [21 Note that 
iteration of the solution here is crucial but p — 2, R — 1 
is approached at late times. 

The numerical results only converge to p = 2, i? = 1 
very slowly. This can also be understood analytically. 
There is a slight correction to this result arising as a 
result of ^ being finite, but it is already insignificant by 
log a ~ 2. A more significant correction to the scaling 
behavior arises as a result of A not being zero. Its size 
can be estimated by assuming that p — 2+e. Substituting 
this into both Eq. l(59|) and Eq. l(60ll we obtain 



'4 



6H^r^^nn,A 



(61) 



Indeed, as f2 grows, the correction decreases until p = 2 
is achieved. In particular, the continued growth of fi 
implies that the anisotropy at the largest scales will con- 
tinue growing as 



(62) 



4- Superhorizon Metric Evolution 



Deriving an evolution equation for the metric high- 
lights the fact that the metric evolution is not sensitive 



to the DGP modification until the acceleration epoch. 
Even in the acceleration epoch the metric obeys a simple 
equation of motion on superhorizon scales. 

By appropriately rearranging the linearized Einstein 
equations and the Bianchi identity, we can obtain a single 
differential equation for the gravitational potentials 



H 



H' 



H'' 



H' 



k 



H' 



H' 
1 



H 



H' 



' 3H2Hr^ 
2Hrr - 1 



9{1 + w)H Hrc 
H' 1 



1 



H 2Hrr 



(w - ci)AE + 



{l + w) 
H' 

""+311 



Se 



(63) 



In the large-scale limit, k/aH 0, Se ~ 0; if we also 
assume that — w or that Ae — 0{k/aH)^<l> from the 
Poisson equation, we obtain 



H' \H' H 



tf = 0. 



(64) 



Note that this equation for the metric perturbations is 
the same as found in general relativity [e.g. [23[ Eq. (50)]. 
Therefore, before the DGP modifications change the ex- 
pansion rate or generate anisotropic stress, the evolu- 
tion of the metric is identical to general relativity. Fur- 
thermore, as pointed out by Bertschinger [2^ in the ab- 
sence of anisotropic stress, ^ = — ^ and this equation 
yields solutions that depend only on the expansion his- 
tory through H even in the acceleration epoch. 

In DGP gravity, however, the anisotropy at the largest 
scales is never negligible and, in fact, grows at late times, 
as discussed in section IIV A 31 We can rewrite Eq. l[64|) 
by defining 



(65) 



In the de Sitter era, when H'/H = and H" / H' = -3, 
and assuming that the two new variables obey a scaling 



solution with ^+ = 
becomes 



and ^_ = A_ 



Eq. ([Ml 



A+aP+ {p\ + 2p+ - 3) + A^aP- {pi + 4p_ + 3) = . 



(66) 



One would expect that would grow with exponent 
p+ = 1, while <P- would decay away with exponent p_ = 
— 1. However, because of the need to preserve the Bianchi 
identity with a non-zero A, as discussed in section [TV A 3| 
the scaling solution is slightly violated. We find that 
p+ = 1 + e while p_ = 1. This leads to the relation 



A- 

a: 



(67) 



where e ^ monotonically. At late times, our solution 
tends to a regime where the ratio 0, with 

<P+ cx a. 



A Limiting Cases and Numerical Solutions 



9 



It is interesting to note that in this opposite Hmit to 
that studied in 



where ^_ < <P+ Eq. also be- 
comes closed and has solutions that depend only on the 
expansion history through H. 




FIG. 4: Ratio of master variable Q for dynamic scaling and 
quasi-static solutions. In QS, f2 responds instantaneously to 
changes in A. Fully dynamic solution to the Bianchi identity 
l|3ip requires time to respond and eventually decays to the QS 
solution. Initializing the calculation at earlier times changes 
neither the scale factor at which f2 responds nor its value 
today. Rapid growth occurs during the time before horizon 
crossing. 



In fact, Eq. l(64|) is equivalent to the statement that the 
Bardeen curvature, ( = <!> + Hq, is conserved during the 
de Sitter era, once the mode leaves the horizon. This is 
also true of the comoving density perturbation, A, which 
saturates to a constant during the de Sitter era. However, 
the comoving and longitudinal hypersurfaces warp with 
a shift Hq that grows without bound. 
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FIG. 5: Evolution of the principle gravitational observable 
$ — ^ {or concordance ACDM, the quasistatic (QS) solution, 
and our new dynamical scaling (DS) solution. For scales 
k > 0.01 Mpc~^ the scaling and QS solutions do not differ 
appreciably: the decay in the potentials is a little faster than 
ACDM as a result of slower growth of density contrast. At 
larger scales, the DS solution exhibits significant additional 
decay owing to the different value of f2 at late times, as ex- 
hibited in Fig.m All potentials normalized to 1 at logo = —2. 




It is useful to summarize the differences between the 
quasi-static (QS) and dynamical scaling (DS) solutions 
uncovered in the previous section. 

Beginning at the initial conditions in the radiation- 
dominated era, the superhorizon value of the master vari- 
able is highly suppressed with H-qs/^qs — O (^{k/aH)^^ 
(see Fig. |4]). As the mode enters the horizon during 
matter domination, the DS solution for J7 grows rapidly 
and then executes damped oscillations around the QS 
solution. This can be understood analytically since the 
Bianchi identity takes the form of a damped oscillator 
in n/a^ that is driven by A. During the time when 
fi significantly deviates from the QS solution, the Weyl 
corrections to the Poisson equation lfT6|) are suppressed, 
since 3> 1, and there is no additional correction to 
the gravitational potentials over and above that of QS 
(see Fig. III). 



• ■ ' 

-2 -1.5 -1 -0.5 

logo 

FIG. 6: Evolution of the comoving density contrast A for 
concordance ACDM, QS and our DS solution. The growth 
function in DGP is suppressed compared to ACDM, even in 
the case of a flat cosmology. There is no significant difference 
between the scaling ansatz and the quasi-static solution, with 
DS departing at most by 2% from QS at the largest scales. 
All quantities are normalized to 1 at log a = —2. 



We find that the results of the scaling solution match 
the quasi-static results for all modes that enter the 
horizon well within the matter-dominated epoch k > 
0.01 Mpc~^ (see Figs. SHl]). For larger scales, this is 
not so: as shown in Fig. lU f2 now only decays toward 
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the QS solution, rather than oscillating around it. Since 
at late times the Weyl perturbations are no longer sup- 
pressed, this now makes a significant contribution to the 
Poisson equation, resulting in additional decay of the po- 
tentials. As shown in Fig. [7l there is a 15% deviation 
from QS in ^ - at the scale k = 0.001 Mpc"^ for the 
chosen sets of cosmological parameters. The direction of 
this effect agrees with estimates made by Lue in j23| and 
is such that the scaling solution is an even worse fit to 
CMB anisotropy data than the quasi-static (see jl^ for 
a discussion). 

On the other hand, up to the present time, the QS so- 
lution for comoving density perturbations. A, is a very 
good approximation for the DS solution at all scales. 
The additional suppression is of the order of 2% for 
k = 0.001 Mpc-i. 
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FIG. 7: Comparison of results of QS and DS solutions: Upper 
panel presents the ratio of # — I?" in the two approximations. 
For modes with k > 0.01 Mpc~^ the two solutions differ by 
less than 2%. The difference is much more pronounced for 
larger scales where fl has not decayed to the QS value, result- 
ing in additional decays of up to 15%. Lower panel presents 
the ratio of comoving density perturbations in the two so- 
lutions: A is affected much less, with approximately a 2% 
deviation from QS at the largest scale, k — 0.001 Mpc~^. 



V. DISCUSSION 

We have introduced a new scaling ansatz which allows 
solutions to Hnear perturbations in the DGP model on all 
scales less than the cross-over scale Tc up to the present 
epoch. The equations of motion for linear perturbations 
on the brane require knowledge of the gradient of the so- 
called master variable into the bulk. The master variable 
obeys a master equation in the bulk. To solve the master 
equation, it is sufficient to have two boundary conditions, 
one on the brane and the other in the bulk. 

Our scaHng solution begins with an ansatz for the 
brane boundary condition: that the evolution of the mas- 



ter variable is scale free on the brane. The second bound- 
ary condition is that the master variable vanishes at the 
causal horizon in the bulk. With these two boundary 
conditions, we solve the master equation to determine 
the gradient. With the gradient known, we can then re- 
place the scale-free ansatz with the dynamical solution 
and iterate the solution until convergence. 

We find that the quasi-static (QS) solution of lUl] is 
rapidly approached once the perturbation crosses the 
horizon. Before horizon crossing there are strong de- 
viations from the quasistatic solution. For modes that 
crossed the horizon only recently during the acceleration 
epoch, we find that the metric perturbation <P — ^ de- 
cays more rapidly that the QS solution. The QS solution 
itself has a stronger decay than the ACDM model. The 
extra decay compared with ACDM is extremely robust 
to changing the gradient of the master variable into the 
bulk, the one variable that is required to close the equa- 
tions of motion on the brane. We consider the observa- 
tional consequences of these results in a companion paper 
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APPENDIX A: SCALING ANSATZ ROBUSTNESS 

Assuming that the scaHng ansatz is appropriate, the 
solution presented in this paper is correct at all scales. 
However, the scale-free assumption for the solution in the 
bulk depends crucially on the existence of only one scale 
in the problem: the Hubble parameter. It is possible 
that through some additional physics in the vicinity of Vc 
this assumption is broken and the evolution of the bulk 
profile will depend on both H and Tc independently. 

However, any violation of scaling will enter the equa- 
tions solely through the off-brane gradient R: this is the 
quantity which we obtained by solving the master equa- 
tion Eq. (|43l) using the scaling ansatz. It is important to 
note that, in our analysis, we have not dropped any terms 
either arising from the Weyl-fluid-driving Bianchi iden- 
tity (|3T|) or the master equation (|27|) . We have assumed 
that the solution in the bulk depends only on yH/£_ but 
the form of the master equation implies that this scaling 
assumption should be a good one. We are only testing 
the robustness of our solution in order to estimate the 
effect of any new physics which might be important at 
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scales around and which is not embodied in the master 
equation already. 

One way of testing the robustness of the scaling ansatz 
is to change the values of R at late times and investigate 
how much of a departure from scaling-ansatz values is 
necessary to significantly change the behavior of observ- 
ables. We concentrate on the change to the evolution of 
the potential — 'F, which drives the ISW effect, as we 
alter the off-brane gradient. Since the QS solution al- 
ready has a significantly sharper decay than ACDM, and 
therefore is a worse fit to the large-angle CMB anisotropy 
[l5| . and the DS solution decays even more rapidly (see 
figure [U , we attempt to violate the scaling solution in 
such a way as to soften this decay. We present the mod- 
ification to R in Fig. [HI we employ a Hnear interpolation 
for R between its scaHng-ansatz value at log a — —2 and 
a chosen off-brane gradient value today, Rq. One should 
note that this breaking is rather extreme, since the scale 
under consideration, k — 0.001 Mpc~^, is inside the hori- 
zon today and it should be well within the quasi-static 
regime at the present time. 



beneath that of the QS solution as Rq is sent to infinity. 

The above considerations show that the new physics 
required around Tc would have to violate the scaling be- 
havior rather strongly in order to give an ISW effect com- 
parable to that of ACDM. 




logo 




FIG. 8: Off-brane gradient for k = 0.001 Mpc"^ for the scal- 
ing solution and the scaling-violating scenarios employed in 
robustness testing. Since any scaling violation is only likely to 
occur at scale factor close to Hrc ~ 1, we modify the gradient 
starting at log a = —2. A value of the gradient today, Ro, is 
chosen and the gradient is interpolated linearly between these 
two scale factors. In this modification, we disregard the fact 
that the mode enters the horizon at late times. 



FIG. 9: Evolution of - •Z' and Z\ for mode k = 0.001 Mpc~^ 
for a selection of scaling-violating scenarios. Increasing Ro 
brings the solution closer to that of QS, and, for very large 
values, reduces the decay of — to that of ACDM. The effect 
of changing Ro on A is much smaller, the quantity remains 
insensitive to the precise details of the scenario. 
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We have found that choosing negative values for Rq 
strengthens the decay of<P — W, with the scaling solution 
approximately replicated for Rq ~ —2. Positive values 
of Ro reduce the rapidity of the decay at late times: the 
QS solution is matched for Rq « 5, which is a value 
much higher than ever achieved by the scaling solution. 
In order to achieve the low levels of decay exhibited by 
ACDM, Rq needs to be set in the vicinity of 50. 

The effect is different for A: here, negative values of Rq 
also decrease the growth rate even further beneath that of 
DS; however, A asymptotes to a value approximately 2% 



FIG. 10: Ratio of <P— tf' and A for scaling and scaling-violating 
scenarios to their values in QS for mode k = 0.001 Mpc~^. 
Dashed line represents the final value of the quantity in QS, 
dotted — in DS, and dot-dashed — in ACDM. The solution is 
very sensitive to negative values of -Ro, but choosing such sce- 
narios only increases the decay, strengthening the ISW effect. 
Positive values of Ro bring the evolution of the observables 
closer to that of the QS solution and, for very large values, 
achieve decays as low as those of ACDM. A is quite insensi- 
tive to the choice of R. All quantities were normalized to the 
same value at log a = — 2. 
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